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Abstract. 

Thermal fluctuations lift the high ground state degeneracy of the classical 
nearest neighbour pyrochlore antiferromagnet, with easy plane anisotropy, giving 
a first order phase transition to a long range ordered state. We show, from 
spin wave analysis and numerical simulation, that even below this transition a 
continuous manifold of states, of dimension JV 2 / 3 exist {N is the number of degrees 
of freedom). As the temperature goes to zero a further 'order by disorder' selection 
is made from this manifold. The pyrochlore antiferromagnet Er2Ti207, is believed 
to have an easy plane anisotropy and is reported to have the same magnetic 
structure. The latter is perhaps surprising, given that the dipole interaction 
lifts the degeneracy of the classical model in favour of a different structure. We 
interpret our results in the light of these facts. 

PACS numbers: 28.20.Cz, 75.10.-b, 75.25.+Z 



In this paper we present a Monte Carlo study and classical spin wave analysis 
of an easy plane antiferromagnet on a pyrochlore lattice. The work was motivated 
by a recent comprehensive study of the rare earth pyrochlore Erbium titanate P^. 
Er 2 Ti 2 07, which orders at 1.2 K [2] has been suggested to approximate the easy plane 
antiferromagnet 0J |S] . We find that the magnetic structure of this simple model 
system is determined by an order by disorder mechanism 6 and that it is the same 
structure as that observed in the experiment. 

The easy plane antiferromagnet 

The pyrochlore lattice (see figure has a rhombohedral primitive unit cell with 
lattice vectors, 



2 2/' V2 2 ' V 2 2 



with a four atom basis at positions, Vi = , V2 = w , V3 = | , V4 = |, as shown in 
figure |2| The units are those of a larger, cubic cell containing four tetrahedral cells 7 . 
The number of primitive cells, N ce u — N/4, where N is the total number of spins and 
the linear dimension of the sample , L = (N ce u/ '4) 1 / 3 . 
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Figure 1: The pyrochlore lattice of corner- Figure 2: Rhombohedral axes with respect 
sharing tetrahedra to the cubic unit cell 



We have studied the D — > — oo limit of the Hamiltonian 

N 



J Si-Sj-D^iSi-Sif 



<i,j> 



where Si corresponds to a local < 111 > anisotropy and the exchange constant J = — 1. 
The spins are thus constrained to a plane oriented perpendicular to the local < 1 1 1 > 
axis: 



S X 1 + Syl 

S X 2 — S y 2 
S x 3 ~ S y 3 
S x 4 ~\~ Sy4 



s z2 = o 
s z3 = o 

S Z 4 = 0. 



(la) 
(16) 
(lc) 
(Id) 



The ground state condition for four antiferromagnetically coupled spins gives 3 further 
constraints [7]: 



S x i + S X 2 + S X 3 + S X 4 = 

Syl + Sy2 + Sy3 + Syi = 

Szl + 5*22 + S Z 3 + S Z 4 = 0. 

In addition the spins are all of unit length: 

ISxl = |S 3 | = |S 3 | = |S 4 | = 1. 



(2 a) 
(26) 
(2c) 

(3) 

Linear equations Qlaf) - ()2cf) can be solved by row-echelon matrix reduction 8 . 
The resulting dependent and independent variables are put into equation © to give 
four equations in terms of five variables. There is therefore one independent and 
continuous degree of freedom for four spins on a single tetrahedron in their ground 
state configuration. Full details of this solution can be found in Ref. [§]• 

The result for one tetrahedron suggests the possibility of a continuous ground 
state degeneracy when the tetrahedra are connected into a pyrochlore lattice and this 
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Figure 3: Msl vs Tj J, showing a first order phase transition to an ordered q=0 state. 





(a) State 1 



(b) State 2 



Figure 4: Different ground states of the easy-plane tetrahedron antiferromagnet. 

is indeed what we found numerically, as shown below. This result is rather different 
from that found in Ref . ^Oj > where the model was first studied. Here the authors only 
identify a discrete degeneracy and it was a surprise to find evidence of fluctuations 
over a continuous set of states. 

Monte Carlo simulations 



Single spin-flip Monte Carlo simulations were performed on systems of size N = 432 
to N — 3456, with periodic boundaries. The simulation lengths were typically 10 5 
moves per spin, with 40, 000 moves used for equilibration. As shown in figure [21 
initial results confirm the strongly first order phase transition presented in Ref. 
The order parameter Msl is the sublattice magnetisation for each of the sites of 
the primitive cell, which we refer to as q — ordering. For T ~ 0.125 J the order 
parameter jumps from a small value characteristic of a finite system in a disordered 
phase to Msl ~ 0.85, indicating an already highly ordered state. 

Closer inspection .1. shows that the ground state chosen by the system is that 
shown in figure 4(a) not that proposed in Ref. jTH] and shown in figure 4(b) The, 



initially proposed state has pairs of spins either aligned anti-parallel to each other or 
perpendicular to each other. We actually find scalar products between spins here of 
—2/3 and 1/3. The spin components of both states are shown in Tabled 

It is easy to show that not all possible ground states have q — ordering. For 
example, in Ref. ^H] 0(L 2 ) states were identified with spins ordered along one axis 
and with disorder in the planes perpendicular to this axis. The transition must 
therefore be driven by an order by disorder [^j mechanism, although no details of 
the entropic forces involved have previously been given. Below the transition it is 
clear that further entropic selection between different q = states must occur as we 



Soft modes in the easy plane pyrochlore antiferromagnet 



4 



Table 1: The spin components of State 1 and State 2. 
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(a) T/J = 0.09 



(b) T/J = 0.05 



(c) T/J = 0.005 



Figure 5: Bond energy distributions P(Eij) for various T/J. 



have already identified two different possibilities in figures 4(a) and |4(b)| while the 
system systematically ends up in state f for T — > 0. 

and 



In figures 5(a) 5(b) 
different q 



5(c) 



we illustrate the evolution of this choice between 
states. We show the distribution of nearest neighbour bond energies 
for three temperatures below the transition, for which, in all cases Msl > 0.9. For the 
highest temperature there is a continuous distribution of values from —I (antiparallel 
spins) to 1/3 where a sharp peak occurs. This figure is consistent with the result 
for a single tetrahedron, that a continuous degeneracy of q — ground states exists. 
One can see that going to the lowest temperatures the distribution slowly becomes 
peaked around energies —2/3 and 1/3, corresponding to the ground state shown in 
figure 4(a) At all temperatures the energy is very close to the ground state energy 



and this figure clearly represents the order by disorder selection of a single point on a 
continuous energy surface. We expect the selection to be due to a singular excitation 
spectrum above the ground state at this point jll) and this is what we show in the 
next section. 



Spin wave calculations 

We consider small fluctuations out of the ground state 1. Working in the laboratory 
reference frame and re-defining Q S, a = 1,4 as the spins in a given primitive cell, we 
make a small rotation, a S x = a S x o + a e compatible with the local planar constraint. 
S x o is the ground state component and a e is the independent variable. a S y ( a e) and 
a S z ( a e) can thus be calculated by application of the constraints [§]■ The Hamiltonian 
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is then expanded to order ( a e) 2 and written in the form 

N ceU »„„ Q=4,/3=4 

H -^ = ^EE E a a< ^ (4) 

i j a — 1.0—1 

where a and /3 denote spins in a unit cell and i and j denote different unit cells 
(since all bonds are counted twice the sum is divided by 2). Note that due to the 
frustration, individual bonds make linear contributions in ( Q e), but as the ground 
state is a global minimum, these sum to zero. Once the dynamical matrix M Q ^(Ry ) 
has been determined, it is Fourier transformed such that : 



to) = w^: E MQ,3 ( R y) ex p (*q- R y), (s) 

and (QJ is reduced to 



M a/3 (_, 

NceU R 



H - H = \ -L- £ E % ^(q) ( 6 ) 

q a/3 

where 

(7) 

V^cell 

M a ^(q) is shown in matrix (JSJl: 

/ 4J 2J(1 + e ( iq ' a) ) \ 

4J 2j(l + e( iq ( c " b ») 

2J(l + e(~ 4q - (c ~ b)) ) 4J 

V 2j(l + e(-"i a )) 4J 

The eigenvalue equation, /(A) = factorizes to the following form: 

(G - A) 2 - iG 2 A) ((G - A) 2 - l -G 2 B^ = (9) 

where G = 4 J, A = 1 + cos(q.a) and 2? = 1 + cos(q.(c — b)), giving 

A(q)±=4j(l±cos(^)) (10) 

X( q ) ± =Aj(l± C o S (^^-)) (11) 



(8) 



One can see that there are planes in reciprocal space perpendicular to a and to 
(c — b) for which A(q) = 0. We have verified that this is not the case for other selected 
q = ground states and it is these soft modes that provide the entropic force for the 
final ground state selection. 

Using the mode counting arguments of Ref. we can determine the effect of 
these zero modes on the specific heat. For a system with no soft modes, each degree 
of freedom will contribute 1/2, in units of ks to the specific heat. In the harmonic 
approximation, the soft modes contribute zero. Assuming that this singularity is 
regularised at higher order by quartic corrections, each soft mode will contribute 
(1/4). 
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0.5 - C t /N 

0.07 T 




2/3 

One plane in reciprocal space contains N ',, zero modes. The number of zero 



modes per unit cell is 2N C 
is thus: 

C h 



-1/3 
ell ' 



since there are two planes. The specific heat per spin 



Nk B T 



C h 



Nk R T 



(4-2AT-f)i + 



2N 



-1/3- 



f f 



8jv 1/3 

JV ceiZ 
1 1 



2 n /3 L 



If i 

11 2 



N 



(12) 

is plotted versus 1/L, a straight line 



where N ce u = 4L d and C is a constant. 

graph should result. This figure has been tested by Monte Carlo (see figure [SJ). Data 
were collected for L = 1 - 7 at a temperature of T/J = 0.0001 with 1, 000, 000 Monte 
Carlo steps per spin and 500, 000 equilibration steps, with five separate simulations 
averaged for each lattice point. There is good agreement with a linear fit, suggesting 
that our order by disorder scenario is correct. However the gradient, while it is the 
correct order of magnitude, is not correct: 0.063 as opposed to the predicted value 
0.0787. We have not, at present been able to explain this discrepancy. Note that 
the entropic contribution to the free energy coming from the soft modes scales as 
N 2 / 3, and is therefore not extensive. This could mean that the ordering within the 
q = manifold occurs at a system size dependent temperature that goes to zero in 
the thermodynamic limit and the discrepancy could come from this effect; although 
the temperature here is very low. More work is required to resolve this point. 



Conclusion 



In this paper we have reported another example of order by disorder in a geometrically 
frustrated magnetic system (see for example ^]^JE]). We have exposed one order 
by disorder mechanism leading to the low temperature selection of an ordered state, 
but there are several interesting statistical mechanics questions that remain open: 
firstly, is there a different, perhaps discrete mechanism, analogous to that in Villain's 
original paper (i . that drives the phase transition, or are the soft modes felt by the 
system, even in the disordered phase? Secondly, are there disordered states with the 
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same soft mode structure? This is the case for the Heisenberg antiferromagnet on the 
kagome lattice ^] EH an d the results here suggest that columnar line defects could 
exist, equivalent to the line defects of the kagome lattice. It would be interesting to 
pursue this question further. 

The project, in this case was strongly motivated by experimental results on 
Er2Ti2C>7 and its real interest lies in the question: is it in fact relevant? There 
are several details that the simple model does not correctly predict; most notably, 
experimentally the transition is second order pQ while here we find a very strongly 
first order transition. Further, experimentally the low temperature specific heat varies 
like T 3 pp, which is indicative of a spin wave spectrum with a quadratic density 
of energetic states. We have calculated the density of eigenvalues above the zero 
frequency modes and find <?(A) = const., which is compatible with having degenerate 
planes in reciprocal space and variations along only one dimension. Here, of course, we 
cannot capture the quantum fluctuations as our model only has one degree of freedom 
per spin. To study this question more closely one would need to add out of plane 
fluctuations, although it is difficult to see how this could lead to a quadratic density 
of states. It might be interesting to look at this problem in the future. However, 
despite these failings our simple model does have one great success: it does predict 
the correct magnetic structure for Er2Ti207 pQ. This seems particularly interesting 
considering that dipolar corrections, which one might expect to be the most substantial 
perturbation ^3] would lead to a different magnetic structure ^5]- A classical model, 
with dipoles would predict the ground state shown in figure |4(b)| which is different 
from the experimentally observed structure. With this result we suggest that our 
model is a relevant starting point to study this compound and it will be interesting to 
see if quantum fluctuations can change the details discussed above. 
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